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Semi-linear response theory determines the absorption coefficient of a driven system using a 
resistor network calculation: Each unperturbed energy level of a particle in a vibrating trap, or of 
an electron in a mesoscopic ring, is regarded as a node (n) of the network; The transition rates (wmn) 
between the nodes are regarded as the elements of a random matrix that describes the network. If the 
size-distribution of the connecting elements is wide (e.g. log-normal-like rather than Gaussian-like) 
the result for the absorption coefficient differs enormously from the conventional Kubo prediction 
of linear response theory. We use a generalized variable range hopping scheme for the analysis. 
In particular we apply this approach to obtain practical approximations for the conductance of 
mesoscopic rings. In this context Mott's picture of diffusion and localization is revisited. 



I. INTRODUCTION 

Semi linear response theory (SLRT) [iHl] provides a 
procedure for the calculation of the absorption coefficient 
of a driven system, assuming that there are well defined 
transition rates Wmn between levels that are ordered 
by energy. In this context it is helpful to regard w^J^ 
as describing resistors that connect nodes of a network, 
a point of view that has become popular in the related 
studies of variable range hopping (VRH) [MS]. In the 
random matrix theory (RMT) framework Wmn is a ran- 
dom matrix whose construction is inspired by analyzing 
the statistical properties of the Hamiltonian of an actual 
physical system lllll^ . Three physical applications have 
been discussed so far: (i) metallic rings driven by elec- 
tromotive force [l5| : (ii) metallic grains driven by low 
frequency radiation Q; (iii) cold atoms that are heated 
up due to the vibrations of a wall [iB] . It is crucial to ob- 
serve that depending on the parameters that define the 
physical model, the matrix Wmn might be banded and 
sparse [Tslfl^ . Consequently, non-trivial results that go 
beyond linear response theory (LRT) are obtained. 

In order to have a precise mathematical definition of 
the RMT model, let us write the random matrix as 



X BiEm-En) 



(1) 



In this expression B{uj) describes the band-profile of the 
matrix, and Xmn is a random matrix whose entries x 
are positive uncorrelated random numbers. If log(a:) is 
widely distributed over many decades, as in the case of 
log-normal or log-box distribution, then we say that the 
matrix is effectively sparse. Sparsity means that the ma- 
jority of elements are very small compared with the av- 
erage value. 

Irrespective of real-space dimensionality, we regard the 
index n of the energy levels as labeling the nodes of a 
ID lattice (see FigH]), hence the Wmn define a ID resis- 
tor network. The inverse resistivity of this network (see 



App.|A| is denoted as w = [[wmn]] and has the meaning 
of diffusion coefficient. In proper units the relation is 

-De = Qe^ X [[W,nn]] (2) 
= G £2 (3) 



where is the mean density of states (DOS). The pa- 
rameter e represents the RMS amplitude of the driving 
field: it is the RMS displacement of a wall element if 
we consider the heating of cold atoms in a trap; it is 
the RMS voltage if we consider a ring that is driven 
by an electro-motive- force (EMF). We assume here that 
Wmn , which holds whenever the standard conditions 
of the Fermi Golden rule (FGR) are satisfied. 

Within the framework of the FGR picture the tran- 
sitions rates Wmn are determined by the matrix ele- 
ments Vmn of the perturbation term in the Hamiltonian. 
The naive expectation is to obtain the Kubo formula 
G = 7rpE((|KrmP))a for the absorption coefficient G. The 
calculation involves a weighted algebraic average 



(IK, 



2^ 



(4) 



where F{uj) is a normalized function that describes the 
spectral content of the driving source, and {■■)uj is defined 
as the average value for (Em— En) uj transitions. A 
more careful inspection reveals that the Kubo calculation 
does not apply to the problem as defined above. In order 
to appreciate the difference we re-write Eq. ^ as 



G 



(5) 



The double average notation indicates a resistor-network 
calculation. This SLRT 'average' is bounded from above 
by the algebraic average of Eq. ^ and from below by 
the corresponding harmonic mean. Later in this paper 
we adopt a generalized variable range hopping (VRH) 
procedure in order to estimate the SLRT average: 



[\Vmn\XF{uj)p- (6) 
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where [..]u: is the typical value for {Em— En) ~ oj transi- 
tions. The notion of typical value will be defined later: 
it is determined by the size-distribution of the matrix 
elements. 

Physically the idea behind SLRT is very simple: in 
order to have "good" absorption it is essential to have 
connected sequences of transitions. Consequently, if Wmn 
is sparse, the traditional Kubo expression provides gross 
over-estimate, because it is based on an algebraic average 
calculation. Consequently, our interest is to calculate the 
ratio between the SLRT and the LRT conductance, which 
we define as the SLRT suppression factor: 



ffSLRT 



((l^™nP». 



(7) 



where ((...))„ denotes the usual weighted algebraic aver- 
age that appears in the Kubo formula. Loosely speaking, 
if the percentage of large in-band elements is s ^ 1, then 
a generalized VRH estimate might lead to a result of the 
type 



const \ 



(8) 



A few publications have been devoted to report various 
partial results that have been obtained using SLRT. The 
purpose of the present paper is to bridge between SLRT 
and the traditional literature, to further develop the an- 
alytical tools, and to provide elaborated tangible results 
that hopefully can be tested in actual experiments. 

Our main focus concerns the Ohmic conductance Gohm 
of small metallic rings, which is related to the G of 
Eq. ([5]) via Gohm = 0eG, where E is the Fermi energy. 
Up to a factor, the perturbation matrix consists of the 
elements Vmn of the velocity operator. Accordingly, Gohm 
is the LRT or the SLRT average over |wmnP- Past liter- 
ature has provided a theory for the conductance in the 
Debye or adiabatic regimes [l7lll8l | where the PGR pic- 
ture does not apply. Diffusive rings have been further 
analyzed [l^l in the Kubo regime, and later weak local- 
ization corrections have been incorporated j20ll2Tj | and 
verified experimentally [23 - |23 |. 

Still neither VRH in real space, nor SLRT response in 
the ballistic regime had been considered in the context 
of mesoscopic conductance. In Pig. [2] we present some 
re-processed numerical results that have been reported 
in Ref . [l^ . These numerical results indicate that indeed 
for both weak and strong disorder the matrix elements 
of the velocity operator become sparse with s ^ 1. As 
explained in Ref. ^15)] this is related to the non-ergodicity 
of the eigenstates. In the present paper we would like to 
present a full analysis of the conductance that starts from 
the strength of the disorder W as an input. The disorder 
determines the sparsity s, and then, using RMT model- 
ing and a generalized VRH approximation, leads to some 
tangible results (Pig. [3|) for the SLRT suppression factor 
5sLRT- We also explain how this factor can be measured 
in an actual laboratory experiment, and how semi-linear 



response can be distinguished from linear response in a 
way that does not involve any ambiguities. 

Outline.— Sec. [ll] motivates the study by introduc- 
ing the physical model, including subsections that re- 
late to the characterization of metallic rings and their 
Kubo-Drude conductance. Some more details are given 
in App. Eland App.O Sec. IIIII discusses the RMT mod- 
eling in general. Sec. llVl brieflv reviews the SLRT calcu- 
lation procedure. Sec. [V] elaborates on the generalized 
VRH approximation. Sec. I VII introduces the analysis of 
some prototype non-Gaussian ensembles. Some more de- 
tails are given in App.lPjand App. [E] Sec. I VIII discusses 
the semiclassical theory of the matrix elements that are 
required for the calculation of the mesoscopic conduc- 
tance. Sec. I Villi discusses the SLRT calculation in the 
ballistic regime. Sec. IIXI discusses the SLRT calculation 
in the Anderson localization regime. Sec. |X] clarifies 
the relation between SLRT and the traditional VRH cal- 
culation. Sec. IXII questions the possibility to get VRH 
from proper LRT analysis. Sec. IXIII contrasts VRH with 
non-thermal hopping due to noisy source. Sec. lXlTIj pro- 
poses how to experimentally test SLRT via conductance 
measurements. Sec. IXIVI summarizes the major obser- 
vations regarding the relation between SLRT, LRT and 
VRH. 



II. PHYSICAL MODEL 

In order to physically motivate the analysis, we con- 
sider a particle of mass m in a rectangular box of length 
Lx = L and width Ly. In one problem, that of Ref.[16|], 
we had assumed Dirichlet boundary conditions and con- 
sidered the response for vibrations of the wall. In the 
present paper we assume ring geometry with periodic 
boundary conditions on L^, and consider the response 
to electro-motive force (EMP). In both cases the Hamil- 
tonian matrix can be written as 



n = diag{^„} + {C/^,„} + f{t){Vrr^,r^} 



(9) 



where n — (n^, Uy ) labels the unperturbed eigenstates of 
a clean box/ring, U{x,y) describes the potential floor 
(either smooth deformation or uncorrelated disorder), 
and V is the perturbation matrix due to the driving. 
Given that the energy of the particle is E we define 
= (2mi;)i/2 and = i2E/my/^. The associated 
number of open modes, i.e. the number of energetically 



allowed ny values, is 



M 



The density of states is 



k^Ly 



m L 

^ 2^^^^^ = ^2^ 



(10) 



(11) 



The static part of the Hamiltonian can be diagonalized, 
and in the new basis the Hamiltonian takes the form 



H = diag{£;„} + fmVrnn} 



(12) 
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where En are the perturbed energies. The power spec- 
trum S{uj) = e'^F{ijj) of the low frequency driving / is ei- 
ther rectangular with sharp cutoff at some frequency Wc, 
or exponential 

n-) = ^oxp(-M) (13) 

We assume that u^c is small compared with any relevant 
semi-classical energy scale, but larger compared with the 
mean level spacing. If the driving is by a thermal source 
then Wc can be identified as the temperature of the source. 
This latter point of view is useful in the discussion of the 
relation between SLRT and VRH. 

In the case of an EMF driven ring e is the RMS of the 
voltage, and the interaction ~](t)V of the particle with 
the magnetic flux f(€) involves V — ~{e/L)v, where v is 
the velocity operator. Hence 

where r is the position operator. Thus an LRT or an 
SLRT study of the conductance reduces to a study of 
the statistical properties of the so called dipole matrix 
elements. These statistical properties become non-trivial 
for either weak or strong disorder, and they should be 
described by a non-Gaussian ensemble. 

The following subsections contains some extra details 
regarding metallic rings, and can be skipped in first read- 
ing. 

A. The characterization of metallic rings 

A metallic ring is characterized by the Fermi veloc- 
ity We, the Fermi momentum ^e, the length of the ring L, 
its width L_L, and the strength of the disorder W. The 
latter determines the mean free path £. The Fermi veloc- 
ity Ue can be regarded as providing conversion between 
"length" and "time" , hence we have two dimensionless 
parameters: the number of open modes A4 ^ {kj^L±)'^^^ 
and the degree of disorder L/£. Formally there is a third 
independent dimensionless parameter k^L, but we as- 
sume it to be very large compared with , and hence it 
has no significant role in the analysis below. The various 
regimes in this problem are described below and in the 
diagram of Fig. 

• Clean ring {L/£) < l/M 

• Ballistic ring {L/l) < 1 

• Diffusive ring {L/£) > 1 

• Anderson regime {L/£) > Ai. 

It is a matter of terminology whether to exclude the 
"clean" case from the ballistic regime, and the "Ander- 
son" case from the "Diffusive" regime. 

The time scale which is associated with the length of 
the rings is t^ = L/ve, the time scale which is associated 



with the scattering is mean free time ti = £/ve, and the 
time scale which is associated with quantum recurrences 
is the Heisenberg time — -Mti,. If the very strong 
condition t^ > t^ is satisfied, then we call it "clean ring" , 
meaning that the disorder does not mix the levels, and 
its effect can be treated using first order perturbation 
theory. 

More generally we define the ballistic regime by the 
condition £ ^ L. If the disorder is strong enough then the 
levels are mixed non-perturbatively leading to genuine 
semi-classical ballistic behavior with 

tL <te < tn [BaUistic] (15) 

In the diffusive regime it is meaningful to define the er- 
godic (Thouless) time via the relation Dot ^ where 
Dq = Ve£, leading to t^^^ — {L/£)tL. In the strict diffu- 
sive regime we have 

te < < in [Diffusive] (16) 

If we have (formally) i^ig > then there is no ergodiza- 
tion but rather a strong (Anderson) localization effect 
shows up. This means that one expects a breaktime ti^^ 
that marks a crossover from diffusion to saturation. A 
standard argumentation (see below) gives the estimate 
tioo = Ai'^ti. One observes that in the Anderson regime 

te < < tn [Anderson] (17) 

The self consistent determination of £^ originates in old 
studies of dynamical localization in the quantum kicked 
rotator problem. Assuming that the localization length 
is £^, the local level spacing is = 7tve/M£^, and hence 
the breaktime is iioc = 27r/Aj. The self consistency con- 
dition is I?<ioc ^ leading to 

£^ w M£ (18) 

The identification of the Anderson regime is via the re- 
quirement L > £^. Finally we note that the largest mean- 
ingful value of disorder is {L/£) = k^L, for which £ equals 
the Fermi wavelength. 

B. The Kubo-Drude conductance 

The coefficient G is defined through the expression 
De — Ge^ for the one-particle diffusion coefficient, where 
e is the RMS of the voltage. Taking Eq. into account 
it follows that 

G-TTfeX (I)' {{\Vrnn\^)) (19) 

The Ohmic conductance is defined as the coefficient in 
the Joule formula GohmC^ for the rate of energy absorp- 
tion. For an Af particle system at temperature T it is 
related to G via a general diffusion-dissipation relation: 

Gohn, = & X G [Fermi] (20) 

Goh.„ = (A/'/T) X G [Boltzmann] (21) 
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The Boltzmann occupation applies to semiconductors, 
where N /L is the density of the particles. For Fermi 
occupation q-e/L is the density of states at the Fermi en- 
ergy per unit length of the ring, and it is in agreement 
with the Boltzmann result if we regard M — g^T as the 
effective number of carriers. 

In the case of a diffusive ring it makes sense to relate 
the diffusion in energy to the diffusion in real space. This 
relation holds in the strict DC limit. Using the Einstein 
relation Gohm = (e/i)^^^^ we deduce that 



V 



TTQe X ((I 

I ))ti'e~0 



(22) 



It is important to keep in mind that for a disconnected 
ring V = but still we can get from Eq. ([2^ a non-zero 
result G ^ Q because the spectral content of the driving 
may have a finite cut-off frequency ujc- 

The reference case for all our calculations is the Drude 
result which is obtained for a diffusive ring in the semi- 
classical approximation (see App.[C]). Assuming a mean 
free path £ we write the Drude result as 



M — 

2nh L 



(23) 



where L is the length of the ring, and Ai is the number 
of open modes (proportional to its cross section) . 

The quantum Kubo calculation gives in leading order 
the same result as Drude: this is well known, and obvi- 
ously it is also a by product of the subsequent analysis. 
One observes that the there is a maximum Kubo conduc- 
tance which is obtained in the limit of a clean ring, i.e. 
for £/L = Ai. For completeness we note that for a ring 
with transmission go, the following formal identification 
applies (see App. [C|): 



L ^ 



90 



90 



<M 



(24) 



This makes transparent the relation between the Drude 
and the Landauer results. 

In later sections our interest is to find the SLRT 
suppression factor (;slrt that determines the ratio 
Gohm/Gnrudc- For this purpose we have to find not only 
the average value of |umnP but also their statistics. 



III. RMT MODELING 

Regarded as a random matrix Vmn is characterized by 
its band profile, and by the size-distribution of its el- 
ements. The standard RMT modeling due to Wigner 
assumes either full or banded matrix with elements that 
are taken out of a Gaussian distribution. But our in- 
terest is in circumstances where the size distribution is 
wide, i.e. the elements of look like realizations of 

a random variable x whose logarithmic value (log(x)) is 
distributed over several decades. 



In practice the IKnnP of a physical model does not 
have an idealized flat band profile. Consequently, we 
write 



IK. 



X C{E,n~En) 



(25) 



where C{lu) describes the band-profile of the matrix. Nu- 
merically the band-profile is obtained by averaging sepa- 
rately each diagonal [{n—m) — const] of the matrix, and 
plotting the result against w — {Em—E„i) ~ {n—m)g^^. 

The question arises, given a matrix Amn that consist 
of real non-negative elements, how to numerically define 
its bandwidth 6, its sparsity s, and the associated dis- 
tribution p{x) of its in-band elements. For the purpose 
of this paper it was important to adopt an unambiguous 
definition of s, which loosely speaking is defined as the 
percentage of large in-band elements. The suggested pro- 
cedure below is based on the participation number (PN) 
concept. The PN of a set {xi} is defined as 



PN = 



(26) 



and reflects the number of the large elements. The pro- 
cedure to determine s and b goes as follows: (1) We 
consider a truncated Amn within the energy window of 
interest; (2) We calculate the band profile by averaging 
separately the elements over each diagonal; (3) We con- 
struct an untextured matrix A^^ by performing random 
permutations of the elements along the diagonals. (4) 
We construct a uniformized matrix A^J by replacing 
each of the elements of a given diagonal by their average. 
(5) We calculate the participation number of the ele- 
ments in Amn ■ This is like counting the number of large 
elements. (6) We calculate the participation number of 
the elements in A™J^ . This is like counting the number 
of in-band elements. (7) The ratio of the numbers that 
have been calculated in the previous step is defined as 
the sparsity s. (8) Likewise the bandwidth b is deduced 
from the number of in-band elements. 

The size distribution p{x) refers to the in-band ele- 
ments. In order to verify that Amn is really like a random 
matrix we perform the SLRT calculation (see Sec. IIV[) 
once on Amn and once on the untextured matrix A^f^ . If 
the results are significantly different we say that texture, 
i.e. the non random arrangement of the large elements, 
is important. The RMT analysis in this paper assumes 
that texture is not too significant. 

In particular we are interested in the bi-modal, log-box 
and log-normal ensembles. The bi-model distribution is 
characterized by the probability p of having a large value 
X — xi, otherwise x — xq ^ xi, hence 



p{x) = (1 - p)S{x — xq) + pS{x — Xi) 



(27) 



In the case of a log-box distribution the variable ln(a:) has 
uniform distribution within [a;o,a;i], hence 



p{x) 



1 



ln(a::i/a;o) x 



(28) 
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In the case of a log-normal distribution the variable ln(x) 
has a Gaussian distribution with mean ln(a;o) and stan- 
dard deviation a, hence 



p{x) 



1 



1 

-e 
2TTa X 



(ln(x/xo))^ 



(29) 



A random variable can be characterized by the algebraic, 
geometric and harmonic averages 



(x) 

exp[(log x) 



(30) 
(31) 
(32) 



The sparsity of a matrix that consists of uncorrelated 
realizations can be characterized by a parameter s or 
optionally by the parameters p and q that are defined as 
follows: 



s 
P 

q 



{xfl{x^) 
Prob(x>(a;)) 

{{x))^^,^J{x) 



(33) 
(34) 
(35) 



By this definition p is the fraction of the elements that 
are larger than the algebraic average and q is the ratio 
between the median and the algebraic average. We re- 
gard a matrix as sparse if s ^ 1 or equivalently if p ^ 1 
or g ^ 1. 



IV. THE SLRT CALCULATION 

As in the standard derivation of the Kubo formula, also 
within the framework of SLRT, the leading mechanism 
for absorption is assumed to be FGR transitions. These 
are proportional to the squared matrix elements |14fmP- 
The power spectrum of fit) is S(u}) = e'^F{uj), where e is 
the RMS value of the driving amplitude. Consequently, 
the FGR transition rates are 



27r- 



{Em— En)' 



;S{E,yi~E„ 



(36) 



From Eqs. ([T]), (1^51) and ([55]) one deduces the identifica- 
tion 



(37) 



The inverse resistivity of the network has the meaning 
of diffusion coefficient, and from the definition of G in 
Eq. ([U we deduce the SLRT formula Eq. Q with 



{{\Vrnn\')) ^ 



F{Em-Er.) 



{Em~E,,y 



(38) 



This should be contrasted with the Kubo formula that 
involves an algebraic instead of SLRT average: 



(39) 



with average over the reference state n. The average is 
done over all the states whose energy En is within the 
energy window of interest. In the metallic context it is 
an average around the Fermi energy. 

It is a simple exercise to verify that if all the 
matrix elements are the same, say iKrinP = cq, then 
((iKrmP)) = Co too. Also it is a simple exercise to verify 
that the SLRT formula coincides with the Kubo formula 
if there is no randomness, i.e. if |KnnP is a well defined 
function of Em— En- But if the matrix is structured or 
sparse then 

mnn\')). < mnn\')) « {{\Vmn\')) . (40) 

If only neighbouring levels are coupled then "adding re- 
sistors in series" (see App. [X| implies equality of the 
SLRT average to the harmonic average: 



{{\Vmn\')).^ 



1 



-\ -1 



Vm 



-F{Em—En) 



(41) 



More generally the harmonic average is a gross under- 
estimate. A generalized VRH scheme that we present 
in Sec. |V] provides the following approximation for the 
SLRT average: 



{{\Vmn\')) 



Wm 



F{u) 



(42) 



where the maximum is calculated with respect to w. The 
typical value | Vmnl^ for w transitions will be defined pre- 
cisely in Sec.|Vl and it reflects the size distribution of the 
matrix elements. The VRH integral Eq. ([6]) is an ad-hoc 
refinement of Eq. ((42|) that better interpolates with the 
LRT result, and therefore it is advantageous for actual 
numerical analysis. 

Further analysis (see Sec. IVip indicates that com- 
pared with the weighted harmonic average ((|KrmP))h, 
of Eq. (j4ip . the corresponding geometric average 



((iKnnI ))g provides in most cases a better lower bound. 



V. THE GENERALIZED VRH 
APPROXIMATION 

A ID network is characterized by its inverse resistivity 
w = [[wmn]]. Inspired by Ref.Q, the inverse resistivity 
can be estimated analytically by finding the maximum 
threshold such that the elements Wmn > wo form a con- 
nected cluster. This leads in the present context to a 
generalized VRH estimate which we explain in the fol- 
lowing paragraph. 

Given a threshold w and truncating the bandwidth 
at w, a sufficient condition for having a connected cluster 
is to have at least one non-zero element per uj segment: 



QeUJ X Prob 



xB{uj)>w larger than unity (43) 



We define the typical value x^j for range uj transitions via 
the relation 



QeUj X FToh{x>Xi^) 



1 



(44) 
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and rewrite the condition for having a connected cluster 
in the foUowing suggestive form: 



(45) 



Thus an under-estimate for the diffusion coefficient is 
D ~ wuP based on hopping rate w with steps uj. The 
VRH estimate is based on the idea to optimize this under- 
estimate with respect to u and w, leading to 



(46) 



where the maximum is with respect to the hopping 
range w. In the FGR context this leads to Eq. (|42|) with 



(47) 



VI. 



SLRT ANALYSIS OF SOME PROTOTYPE 
NON-GAUSSIAN ENSEMBLES 



In this section we derive results for SLRT suppression 
factor gsLRT for the bi-modal, for the log-box, and for 
the log-normal ensembles. The bi-modal distribution is 
the simplest for pedagogical purpose, while the log-box 
and log-normal ensembles are of greater physical rele- 
vance. The main results are summarized below, while 
further details of the calculation are given in App. iPjand 
App.|El Fig.|3]presents the outcome of numerical analysis 
that tests the accuracy of the generalized VRH approxi- 
mation. In later sections we shall see that the presented 
results are of relevance to the study of conductance in 
the limits of strong and weak disorder. 

The bimodal ensemble.— In this case there is a 
minority of large elements {x — xi) that have percentage 
p <C 1, and a majority of small elements (cc = <C 1) 
that have percentage 1— p. Consequently, the typical 
value for w transition has a percolation-like crossover 
from Xuj — xq to x^j — xi at the frequency lo = (PeP)^^- 
Therefore, 



gF(a;-0), bp<l 
[l/p) F{1/{q^p)), bp>l 



(48) 



where b — g^Uc is the dimcnsionlcss bandwidth, and 
q ~ xo/{pxi). The first expression reflects the possibil- 
ity of majority dominance of the small elements, while 
the second expression reflects the possibility of minority 
dominance of the large elements. Note that the VRH 
approximately implicitly assumes that F{ljj) drops (say) 
exponentially such that gsLnx ^ 1, otherwise the result 
cannot be trusted. 

The log-box ensemble.— In this case the probabil- 
ity distribution of log(x) is uniform over many decades. 
Therefore, it is reasonable to assume that the result for 
5sLRT is minority dominated. It is natural to character- 
ize the log-box distribution of Eq. (|28)) by a parame- 
ter p = (ln(a;i/a;o)) ^ ^ and to realize that the percentage 



of large elements is p « —pin p. Note that the corre- 
sponding sparsity parameter is s « 2p. The typical value 
for UJ transitions is 



1 

- exp 

P 



1 



PQ-eUJ 



((^))a 



(49) 



and the VRH estimate, assuming an exponential band- 
profile gives 



<7SLRT 



1 

- exp 

P 



1 

p b 



1/2- 



(50) 



Note the similarity, as well as the subtle difference, com- 
pared with the bimodal minority dominance expectation. 

The log-normal ensemble.— In this case the proba- 
bility distribution of log(a;) is a Gaussian centered around 
the median. Therefore, it is reasonable to assume that 
the result for (^slrt is majority dominated. It is natural to 
characterize the log-normal distribution of Eq. (|29l) by a 
parameter g, which is defined as the ratio of the median 
to the algebraic average. Note that the corresponding 



sparsity parameter is s 
the result 



The VRH calculation gives 



5SL 



q exp 



factor X In ( - ) In (6) 



1/2 



(51) 



where the factor is determined by the bandprofile (it 
is 2 for an exponential bandprofile and 4 for a rectan- 
gular bandprofile). Note that ^slrt ~ Q is the simplest 
guess that reflects the majority dominance expectation. 



VII. THE SEMICLASSICAL ESTIMATE OF 
THE OHMIC CONDUCTANCE 

There is a well established semi-classical procedure 
to deduce the algebraic average of the matrix ele- 
ments IKnnP that correspond to the energy difference 
Lo = Em— En from the associated correlation function 
(y{t)V{Q)). We would like to apply this procedure in or- 
der to estimate the conductance of metallic rings. Hence 
our interest is in the matrix elements of the velocity op- 
erator. The semi-classical estimate is based on the fol- 
lowing observation: 



-FT 



{v{t)v(Q)) 



(52) 



where FT stands for Fourier transform. The velocity- 
velocity correlation function can be obtained via a 
time-derivative of the time dependent diffusion coeffi- 
cient I'(i), which is the time derivative of the spread- 
ing((r(0-r(0))2). 

In the Drude "classical" approximation one assumes 
an exponential decay of the velocity-velocity correlation 
function, and long time diffusion I?o as determined by the 
mean free path (see App. [C]). This is satisfactory in the 
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ballistic and diffusive regimes, and leads to a Lorentzian 
line shape: 



{(KM. 



1 



Riuj) 



b ""l + iteujy 



(53) 



where b = A4L/£ is the dimensionless bandwidth of the 
matrix, and R{uj) = 1. 

But in the Anderson regime we know that there is a 
breaktime ^lo^ that marks the crossover form diffusion to 
localization, and hence for a bulk system formally V — 0. 
Consequently, the FT consideration of Eq. (|53p leads to 
the conclusion that in the limit w — >■ the band-profile 
should vanish if the system is infinite. The simplest rea- 
soning leads to the expression 



i?(w) 



1 



1 + (iiocw)- 



(54) 



where the first term reflects the finite length of the system 
and it is deduced using Eq. We shaU see in Sec lIXl 

using a different more refined approach, that the second 
term is almost correct. Namely, the more careful analysis 
using the Mott's picture predicts that the small frequency 
dependence is not [w]^ but [wlog(w)]^. 

The quantum mechanical analysis should further take 
into account (i) the statistics of the levels and (ii) the 
fiuctuations in the size of the matrix elements. The for- 
mer implies wiggles in i?(aj) for small frequencies, while 
the latter imply that the average size of the matrix el- 
ements does not necessarily refiect their typical value. 
Fig. 2] summarizes the dependance of the matrix elements 
on the disorder. 

It should be clear that we always have the sum rule 



El 



(55) 



In the clean ring limit the sum is dominated by the di- 
agonal or near diagonal element, while all the other off- 
diagonal elements become negligible. Still the estimate 
Eq. ([5^ for the other off-diagonal matrix elements re- 
mains valid and can be justified using 1st order pertur- 
bation theory. If the ring is ballistic (but not "clean") 
then the semi-classical estimate Eq. ([55)1 implies that the 
large elements form a band of width 6 > 1. If the matrix 
is not sparse, then the contribution of all the b in-band el- 
ements to the sum rule is comparable. But if (say) only 
a fraction s <C 1 of of elements are contributing, then 
their typical value |wmnP ~ {{\vmn\'^)) / s is much larger 
compared with the average. We shall come back to a 
more detailed discussion of 'sparsity' in the subsequent 
sections. 

In the diffusion regime R{uj) mainly reflects the level 
spacing statistics of the individual levels, which is a "mi- 
croscopic" effect that leads to small weak localization cor- 
rections that had been studied extensively (20ll2l| . But in 
the strong localization Anderson regime the implication 
of the breaktime leads to the dramatic conclusion that 



R{u!) ^1 for w <C A^, where the local level spacing 
is not related to the volume dependent microscopic level 
spacing g~^, but to the strength of the disorder. 

In the Anderson regime it is evident that ((|wmn P))tj 
is not the typical value of the matrix elements. Roughly 
speaking, and disregarding the uj dependence. 



(56) 



TT exp --r- 



M 



If 



where r € [0, L/2] has a uniform distribution, implying 
a log-box distribution for the size of the elements. Ac- 
cordingly, the typical value is exponentially small in the 
length of the ring, while the average is determined by 
the small percentage of large elements, and comes out in 
agreement with the semi-classical estimate. In Sec. lIXI we 
further elaborate on the statistical analysis of the spar- 
sity in the Anderson regime using the Mott's picture of 
localization. 



VIII. THE RMT STATISTICS IN THE 
BALLISTIC REGIME 

For zero disorder W = Q each energy level is doubly 
degenerate in the basis of real eigenfunctions, and the 
couplings are pairwise, i.e. the matrix element between 
states of different energies is zero. See Fig. [5l Conse- 
quently, the EMF cannot induce connected sequences of 
transitions, and the SLRT conductance should be zero. 
The non-zero elements of the perturbation matrix ac- 
cording to the sum rule (Eq. [55]) are |u,i„i| = v-^. The 
algebraic average of the near diagonal elements equals 
this value (of the large size elements) multiplied by their 
percentage po ~ 1/2. Consequently, 



(57) 



For sufficiently small W these large size matrix elements 
arc not affected, and therefore, the algebraic average 
stays the same. Consequently, in the clean ring limit 
the Kubo conductance is formally finite, and attains the 
maximal value as discussed with regard to Eq. ([221) • 

In the clean as well as in the whole ballistic regime the 
algebraic average {{\vmn'\^))ui does not reflect the sparsity 
and the textures of the Vmn- See Figs. [S] and [T] When 
we look on the image of Vmn the immediate reaction is 
to be impressed by the texture, and therefore we discuss 
it flrst. Subsequently, we discuss the sparsity, which is in 
fact more signiflcant for the analysis. 

The mean DOS of the 2D ring is Qs^. But ^ Ly 
and, therefore, it is not uniform. As the disorder W is 
increased, levels start to mix flrst in the high DOS re- 
gions, and only later in the low DOS regions. This is 
the reason for the appearance of textures. Let us be 
more detailed about the non- uniformity of the DOS. As 
a function of the energy E each time that a mode is 
opened the DOS is boosted. Consequently, is modu- 
lated. This systematic modulation is associated with the 
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opening of a single additional mode at every threshold 
energy and, therefore, scales like l/M.. On top there is 
an additional weaker non-systematic modulation of the 
DOS, because the levels of low density modes add up to 
the levels of the high density modes. It is the latter type 
of modulation which is reflected in Figs. [5] and [U where 
the energy window contains throughout exactly 10 open 
modes. 

In the regions where levels are not yet mixed one can 
estimate the majority of small matrix elements using first 
order perturbation theory: Due to the first-order mixing 
of the levels, the typical overlap |(m|n)| between per- 
turbed and unperturbed states is 



\{m\n) 



En — En 



(58) 



The typical size of a small Vnm element is the multipli- 
cation of this overlap, calculated for {En—Em) ~ Q^~^, 
by the size of the non-zero \vnm\ = Ve element. As a- 
priori expected this first order estimate gives a result that 
agrees with the semi-classical estimate Eq. ([55)1 evaluated 
for Lo g~^. Thus 



[for white disorder] 



(59) 



Above some threshold, first order perturbation theory 
fails everywhere, meaning that non-perturbative mixing 
takes place in any energy. Still, due to the modula- 
tion of the DOS, the mixing range is wider in the near- 
thresholds energies, and therefore the matrix elements 
there are smaller. So now we have the opposite situation, 
of high DOS bottleneck instead of low DOS bottleneck. 

One easily observes that the crossover from weak 
disorder (that features separated mixing regions and 
low DOS bottlenecks) to stronger disorder (that features 
a connected mixing region and high DOS bottlenecks) 
is associated with the crossover from the "clean" to the 
"ballistic" regime. The width of the crossover region de- 
pends on the non-uniformity of the DOS, and therefore 
diminishes as the number of open modes becomes large. 

The above reasoning implies that the texture might be 
important in the SLRT analysis primarily in the clean 
ring regime, but much less in the genuine ballistic regime. 
But what about sparsity? Using the FGR in order to de- 
termine the energy range over which mixing takes place, 
we obtain an estimate for the bandwidth of the pertur- 
bation matrix 



b = 2n gl\U„ 



(60) 



which agrees with the semi-classical estimate. But in the 
ballistic regime b < M. This means that a typical eigen- 
states cannot occupy all the A4 open modes. Rather it 
has there a participation number M = b smaller than M 
(see Fig. [2]). Consequently, we deduce that the sparsity 
of the perturbation matrix is 



M 
M 



L 



where the identification of s with q"^ is based on the as- 
sumption of a log-normal distribution which we further 
discuss in the next paragraph. Unlike the texture, the 
sparsity persists via the whole ballistic regime up to the 
border with the diffusive regime. For this reason we re- 
gard the sparsity as the main ingredient in the SLRT 
analysis. 

The discussion of sparsity in the previous paragraph 
is somewhat meaningless unless one specifies the distri- 
bution to which s refers. At this point of the discussion 
it is essential to distinguish between white disorder for 
which the scattering is isotropic, and smooth disorder 
for which only nearby modes are coupled (small scat- 
tering angle). The latter applies if the potential floor 
within the ring has a smooth rather than erratic variation 
with respect to the Fermi wavelength. Assuming smooth 
disorder it becomes essential to extend the perturbation 
theory of App. [B] beyond first order. It makes sense to 
say that |wm„| ~ \W\^, where the order r is bounded 
by M. Therefore, log(|wm„|) has some bounded distri- 
bution which can be approximated (say) by a Gaussian. 
It follows that a log-normal ensemble should be qualita- 
tively appropriate to describe the statistical properties. 
It follows that (7sLRT can be estimated using Eq. ([STjl with 
the q of Eq. ([61]) . On the other hand in the case of 
white disorder \vmn\ of the majority elements is given by 
first order perturbation theory, and then one should use 
Eq. §Q with the q of Eq. 



IX. 



THE RMT STATISTICS IN THE 
ANDERSON REGIME 



The simplest picture of localization regards the lat- 
tice as composed of segments of size £{, and assumes 
that each eigenstate is well localized in one of this seg- 
ments. Accordingly, non-negligible matrix elements are 
only between states that reside in the same space seg- 
ment. We shall refer to this as the "zero order" picture. 
Taking into account that the matrix element of the ve- 
locity operator are related to those of the position op- 
erator by the relation = i^^kmnP, it follows that 
R{uj) ^ (w/A^)^ in consistency with the semi-classical 
reasoning of Sec. IVIIl which is summarized by Fig. |8l 

In order to refine this picture we use the following pro- 
cedure due to Mott. The zero order basis is determined 
by ignoring the possibility of the particle to hop from seg- 
ment to segment. In order to find the "true" eigenstates 
we have to take into account the residual interaction. It is 
reasonable to postulate that if the distance between two 
zero order eigenstates is r = r„ — 
interaction is 



then the residual 



= Aecxp(-|r|/£e) 



(62) 



[for smooth disorder] (61) 



The prefactor is the natural educated guess, which is later 
justified (see below) by requiring consistency with the 
semi-classical result. 
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If we have two zero-order eigenstates that do not re- 
side at the same segment, but have distance r in space 
and distance e in energy, then the true eigenstates have 
energy difference lo — -f k^, and the dipole ma- 
trix element becomes \r„in\ = (k/uj) x (r/2) instead of 
zero. OriginaUy the zero-order eigenstates had a den- 
sity {qe/ L)dedr, but now the region |e| < exp(— |r|/^^) 
is depleted, and forms a density dcu/ of so called Mott 
resonant states. If we slice all those states that have en- 
ergy difference uj then 



where 



r e-l' 



off res. 
on res. 



(63) 



(64) 



This implies that the size distribution of the \vmn\ ele- 
ments that reside inside a band of width ui is within 



M 
M 



L 



1 



for |w| > (65) 
for |w| < A^ (66) 



Compared with Eq. (154)) this is a refinement that takes 
properly into account the w dependence of the matrix 
elements. Disregarding a logarithmic correction it repro- 
duces the semi-classical result Eq. ([55]) . 

If we ignore the Mott resonant states, then a log-box 
distribution is implied. The Mott resonant states form 
a box distribution on top. In a log-scale the Mott reso- 
nant states contributes a peak of large elements. But this 
peak does not affect the x^ calculation. Consequently, for 
practical purpose we can regard the matrix elements in 
the SLRT calculation as having a simple log-box distribu- 
tion as reflected by the crude approximation of Eq. (1561) . 
The sparsity of this distribution is characterized by 



Mi 



(67) 



and the SLRT suppression factor is given by Eq. (|50|) . 



X. SLRT VS VRH CALCULATION 

In order to appreciate the similarities and the dif- 
ferences between SLRT and the conventional Hopping 
calculation, we cast the latter into the SLRT language. 
Eq. (4.4) of Ref. 7] for the DC Hopping conductance due 
to phonon induced transitions is 



to the Fermi occupation factor, the network contains ef- 
fectively Af — QeT nodes. The division by Af is required 
because we have defined the [[...]] as inverse-resistivity 
and not as inverse-resistance of the network. 

The occupation factor {\-f{En))f{E^)/T gives 0{\) 
weight only to the Af levels that reside within a win- 
dow of width T. If we ignore the relaxation effects 
and regard the fluctuating environment as a noise source 
that induces transitions oc expdE'm— -E'nI/r), we 

still should get the same result for G, even if we omit the 
occupation factor. This point of view allows to bridge be- 
tween the noisy driving problem that we consider in this 
paper and the phonon-induced hopping in the prevailing 
literature. 

The Einstein relation Gohm = / L)^ q-i^T) relates the 
conductance and diffusion in real space. We deduce that 



(69) 



This should be compared with the SLRT expression for 
the noise induced energy diffusion 



= 



11 , 

-^mnJ J _L 



(70) 



Here the resistance of the network is calculated between 
states that reside far away in energy. The SLRT result 
for I?E and the hopping implied result for 2? are both 
simple and manifestly equivalent: The diffusion coeffi- 



cient equals the transition rate 



J] times the step 



squared. In the SLRT calculation the step in energy 
space is I/^e, while in the standard real space analy- 
sis the step is L/Af. Optimization of the hopping with 
respect to the distance w in energy is equivalent to opti- 
mization with respect to the distance r in space. 



XI. CAN WE GET VRH FROM KUBO? 

The analysis that we have introduced in this paper 
gives the impression that SLRT is essential in order to 
derive the VRH result. This statement looks to be in 
contradiction with the prevailing common wisdom, and 
therefore deserves further clarification. In the discussion 
below we explain that VRH can be obtained from Kubo 
for an artificial toy model, but not for the physical model 
that we have analyzed in this paper following Anderson 
and Mott. 

It is instructive to point out that the Kubo formula 
Eq. ([5^ can be rephrased as saying that 



Gohn 



I 

77 



T 



(68) 



The notation [[...]1|| implies that the resistance of the 
network is calculated between states at the same energy 
E ^ Ee, that reside in opposite sides of the sample. Due 



D 



(71) 



where D„ is the diffusion coefficient for a spread- 
ing process that start at state n. If we consider an 
artificial model where the eigenstates are labeled as 
n = (iv), with energies Eii, = and matrix elements 
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exp( 



) 1/^4)7 such that 



exp 



4 



(72) 



then all the Dn are the same value. Furthermore, their 
common value is given by a VRH-like expression which 
reflects an optimization of the hopping distance. Conse- 
quently, the average D is also given by the exactly the 
same VRH-like expression. 

However, in the physical model that we have consid- 
ered in this paper the Z?„ in the Anderson regime are typ- 
ically dominated by one term only, and therefore wildly 
fluctuate. It is then clear that an algebraic average would 
give a very large result which is dominated by the minor- 
ity of large elements. In fact our analysis, which merely 
reproduces Mott's original analysis, shows that up to log- 
arithmic correction the Kubo formula gives G oc w^. In 
order to get VRH we have to perform an SLRT analysis 
rather than LRT analysis. 

In the above discussion one could wonder whether a 
good strategy for obtaining an SLRT estimate would be 
to take a harmonic instead of algebraic average over Z?„ . 
In fact there are circumstances where such procedure 
gives a very good result [l[. However, in general such 
procedure is expected to underestimate the correct re- 
sult, because it is based on the assumption that the hop- 
ping is always with the same optimal step, as in series 
addition of resistors, without the possibility to bypass in 
parallel. 



XII. VRH VS HOPPING 



These can be adjusted such that FGR transitions are the 
dominant mechanism for energy absorption. This ex- 
cludes the adiabatic regime where near nei ghb or transi- 
tions dominate either due to Landau-Zener [l8| or Debye 
relaxation mechanism |17i |. 

Assuming that FGR transitions are the dominant 
mechanism, this does not automatically imply linear re- 
sponse. The rate of the driven transitions can be smaller 
or larger compared with the environmental induced rate 
of transitions, and accordingly we expect a crossover from 
LRT to SLRT [1]. 

The simple minded indication for semi-linear response 
is a drop in the value of the absorption coefficient if the 
driving is strong enough (see estimates below) . What can 
be measured is the SLRT suppression factor (7slrt and its 
dependence on the spectral content of the driving. 

As observed in Ref.Q, the distinction of semi-linear 
from linear response is not ambiguous. The theory is 
called SLRT because on the one hand the power spectrum 
S{lli) I— \S{ijj) leads to Z) \D, but on the other hand 
S{uj) H> Si{uj) + S2{uj) does not lead to D ^ Di+ D2. 
This semi-linearity can be tested in an experiment in or- 
der to distinguish it from linear response. 

Let us discuss in more details the experimental con- 
ditions that are required in order to observe semi-linear 
response. The problem is characterized by the following 
parameters: 



system 
driving 
bath 



(wc,£) 

(70,7rlx) 



(75) 
(76) 
(77) 



It is customary to assume that a noisy non-thermal 
source has a Lorentzian power spectrum: 



1 



F{u,) = - 



(73) 



Let us consider the Anderson regime and assume that 
Wc <C Aj. It should be clear that the VRH result is not 
applicable here. This is because the transport is domi- 
nated by w > transitions. In this case SLRT give the 
same result as Kubo, which we call simple hopping [To| : 



V 



(74) 



where tc = '^j^c- This is as expected from heuristic con- 
siderations. It describes a random walk hopping process 
with steps of size £^ and time . This type of result has 
been highlighted in old studies of the quantum kicked 
rotator problem (STj . 



XIII. EXPERIMENTAL DEMONSTRATION OF 
SEMI LINEAR RESPONSE 

For a given metallic ring the experimentalist has con- 
trol over the frequency and on the strength of the driving. 



where 0;°^" is the frequency that characterizes the semi- 
classical motion; ujq = g^^ is the frequency corresponding 
to the mean level spacing; Wc and s are the cut-off fre- 
quency and the RMS value of the driving (EMF); and 
70, 7rix are the dephasing and the relaxation rates due to 
the environment. 

As already stated we are not interested in adiabatic 
driving (wc < ujq) but rather in what we call DC driv- 
ing. The conditions that have to be satisfied in an SLRT 
oriented experiment are: 



DC driving 
FGR condition 
LRT condition 
SLRT condition 



Wo < < ljI"" 

OJQ OJc 
We ^ 7 
7 < We 



(78) 
(79) 
(80) 
(81) 



where is the FGR transition rate (Eq. (1551) ). 

There are several experimental methods which could 
support the theoretical predictions of our paper. The 
experiment can be based on metallic rings (gold 25lf26j . 
copper [28j, silver [29ll ) . GaAs and other semiconductor 
heterostructures (|27||30|). molecular wires, etc. To esti- 
mate the experimental numbers let us consider a semicon- 
ductor (GaAs) ring driven by time-dependent magnetic 
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flux 



M = 5, L = 0.1 fim, £=50/iTO (82) 
= 2.7 X 10^ m/s (83) 

The long mean free path is required in order to be deep 
in the ballistic regime with sparsity 



M- 



0.01 



(84) 



By Eq. (|lip the mean level spacing is 
2v^ 



Wo 



ML 



1 meV 



(85) 



The ballistic time is — L/vp « 3.7 x 10 s hence 
27rwp 



L 



11 meV 



(86) 



which is ^ 10^'^ i/z in frequency units. In order to satisfy 
the DC driving condition we assume a power spectrum 
of width Wc < ujl^'- The EMF is induced by a time- 
dependent magnetic field e sa [ijj^L?)B. The FGR rate is 
estimated using Eq. p6p with En-^Em ^ ujq 



We 



-2 _ c ^ ^2 



X £ 



(87) 



In order to satisfy the FGR condition the magnetic field 
should be at least 180 G. The expected crossover be- 
tween linear to semi-linear response occurs for ~ 7. 
Assuming for example t^^^ ~ 2 x 10~^^ s we get 7 ^ Sujq 
leading to 



320 G 



(88) 



Under the above conditions we expect that as B is in- 
creased there will be crossover from linear to semi-linear 
response with suppression factor (7slrt ~ 0.3, where we 
used Eq. (|5T|) . The crossover is of course not sharp be- 
cause We is after all distributed over a wide range. In fact 
the functional shape of the crossover can be used in order 
to deduce this distribution ^38] ■ In any case it should be 
re-emphasized that the experimental verification for the 
nature of the crossover requires merely to test whether 
the absorption rate depends in a non-linear way on the 
spectral content of the driving. 



XIV. SUMMARY AND DISCUSSION 

Possibly the nicest thing about SLRT is that it consists 
a natural extension of LRT, that places under one roof 
various results for the conductance in different regimes. 
It should be clear that in the strict DC limit (wc ^ 0), ir- 
respective of the functional form of the power spectrum, 
we always get for G a result that formally agrees with 
the Landauer formula. See the discussion in Sec. Ill Bl In 
the diffusive regime it becomes equivalent to the Drude 



formula with small weak localization corrections. But in 
the other regimes (Anderson, Ballistic), if the low fre- 
quency driving has some arbitrary spectral content, then 
very different results are obtained (Hopping, VRH, gen- 
eralized VRH). 

It is interesting that in our "minimal" treatment of the 
problem there is no need to introduce relaxation due to 
phonons in order to get a VRH result. Rather, we regard 
VRH as arising from the competition between the sta- 
tistical properties of the matrix elements and the power 
spectrum of a noisy driving field. 

The formalism allows to take various limits involving 
the size of the system (L), the driving frequency (wc) 
and its intensity (e), and the rate of the environmentally 
induced transitions (7). The order of the limits is very 
important. In particular: if we take the limit L — 00 
followed by £ ^ 0, keeping 7 constant, then we get LRT; 
while if we take L — > 00 followed by 7 — 0, keeping e 
constant, then we get SLRT. Also note that if we keep L 
constant and take Wc — >■ we get the adiabatic limit and 
not the DC limit of LRT / SLRT. 



We have dedicated Sec. IXIIII to introduce actual es- 
timates that are required in order to observe SLRT in 
a real experiment. It is important to realize that the 
experimental procedure allows to distinguish in a non- 
ambiguous way between LRT and SLRT by playing with 
the spectral content of the driving source. Furthermore, 
one can test specific predictions for the ^slrt suppression 
factor, e.g. Eq. dSD]) with and Eq. ^ with ^ or 
(|FT|) . We note that the explicit incorporation of the envi- 
ronmentally induced transitions into the resistor-network 
calculation, and the subsequent analysis of the resulting 
SLRT steady state is straightforward [38| . 

The SLRT calculation is based on a resistor network 
picture of transitions between energy levels, for which 
an RMT framework is very appropriate and effective. In 
the so called "quantum chaos" context Wigner (in the nu- 
clear context) and later Bohigas (in the mesoscopic con- 
text) have motivated the interest in Gaussian ensembles, 
but there are circumstances where non-Gaussian ensem- 
bles are appropriate, which lead to novel physics. Indeed 
we have faced in this paper the analysis of log- normal and 
log-box ensembles corresponding to the weak and strong 
disorder limits. We have demonstrated that for such en- 
sembles a large SLRT suppression effect is expected, that 
could not be anticipated within the LRT framework. 
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Appendix A: The resistor network calculation 

In this appendix we explain how the inverse resis- 
tivity G = [[GnmW of a one-dimensional resistor net- 
work is calculated. We use the language of electrical 
engineering for this purpose. In general this relation is 
semi-liner rather than linear, namely [[AG]] — X[[G]], but 
[[A + B]] ^ [[A]] + [[B]]. The experimental implications 
of this observation in the SLRT context are discussed in 

SecKml 

There are a few cases where an analytical expression is 
available. If only near neighbor nodes are connected, al- 
lowing Gn,n+i = g-n to be different from each other, then 
"addition in series" implies that the inverse resistivity 
calculated for a chain of length N is 



G 



1 ^ 1 

-T- 



(Al) 



If Gnm — gn-rn IS a function of the distance between the 
nodes n and m then it is a nice exercise to prove that 
"addition in parallel" implies 



G 



M <^ L, and lattice constant a. The longitudinal and the 
transverse hopping amplitudes per unit time are cy and 
c_L , respectively. The random on-site potential in the An- 
derson tight binding model is given by a box distribution 
of width determined by W. 

The numerical calculations of Ref.fisl] assume C|| = 1 
and c_L = 0.9. Thus in the middle of the band there is a 
finite energy window with exactly A4 = AI open modes. 
Rings of length L = 500 with M = 10 modes has been 
considered. In our re- pro cessed Fig. [2] the default cutoff 
is QeCiJc ~ 7 as in Ref.[l5|, but as the disorder becomes 
weaker it is adjusted such that the DC condition Pe^^c ^ b 
is always satisfied. 

For white (uncorrelated) disorder the Hamiltonian is 
given by Eq. ([9]) with the isotropic scattering term 



\Ur, 



ML 



(Bl) 



The eigenstates of the Hamiltonian can be found nu- 
merically. The degree of ergodicity is characterized by 
the participation number FN = [XI P^]"^; which is cal- 
culated in various representations: in position space 



Pr 



— \{fx,ry\'^)\'^ , in position-mode space pr^ 



\{rx, fcjyl*)!^, and in mode space pk^ = I]^^ {{r^, ky\^)\'^ , 



(A2) where ky = [7r/(Al+l)] x integer. 



In general an analytical formula for G is not available, 
and we have to apply a numerical procedure. For this 
purpose we imagine that each node n is connected to a 
current source J„. The Kirchhoff equations for the volt- 
ages are 



(A3) 



This set of equation can be written in a matrix form: 

GV = I (A4) 

where the so-called discrete Laplacian matrix of the net- 
work is defined as 



Gn 



(A5) 



This matrix has an eigenvalue zero which is associated 
with a uniform voltage eigenvector. Therefore, it has 
a pseudo-inverse rather than an inverse, and the Kirch- 
hoff equation has a solution if and only if X^n ^" ~ 0. In 
order to find the resistance between nodes n^^ — and 
J^our = N , we set /q = 1 and /at = — 1 and /„ = other- 
wise, and solve for Vq and Vn- The inverse resistivity is 
G=[{V^^Vn)/N]-\ 



Appendix B: Model details 

In the numerical study we consider the Anderson tight 
binding model, where the lattice is of size L x M with 



Appendix C: The Drude formula 

The velocity-velocity correlation function, assuming 
isotropic scattering, is proportional to the survival proba- 
bility P{t) = e~*/** . Ignoring a factor that has to do with 
the dimensionality c? = 2, 3 of the sample the relation is 



{v{t)vm « vi Pit) 



(CI) 



The rate of the scattering can be calculated from 
the FGR, also know as the Born approximation 



^ = 27r£)E I C^mn I 



= —W 



(C2) 



where in the last equality we used Eq. (jBl[) . From here we 
deduce that the mean free path (disregarding prefactors 
of order unity) 



a \w) 

and the diffusion coefhcient in real space 

1 
2 



Vo = 



{v{t)v{0)) 



(C3) 



(C4) 



By the Einstein relation we deduce the Drude formula 

, 2 ^2 



Gc 



(i) 



A literally equivalent route to derive the Drude formula 
is to semi-classically deduce {{\vmn\'^}} from the velocity- 
velocity correlation function as in Sec lVIIl and then to 
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substitute in Eq. (IT9|) . This has the advantage of aUowing 
easy generahzations of the Drude formula in the baUistic 
and in the Anderson regimes. In this context it is useful 
to realize Q that in the semi-classical picture the integral 
over the velocity-velocity correlation function is related 
to the transmission go of the ring (if it were dissected). 
This leads to the identification in Eq. ([M)) . 

In the diffusive regime Mott has demonstrated that 
it is optionally possible to obtain a direct estimate of 
the dipole matrix elements, using a random-wave pic- 
ture. Namely, it is assumed that in the diffusive regime 
the eigenstates of the Hamiltonian are ergodic in posi- 
tion space, and look like random waves with a correlation 
scale £. Locally the eigenstates are similar to free waves. 
The total volume L'^ is divided into domains of size C^. 
Hence we have (L/i)'^ such domains. Given a domain, 
the condition to have non- vanishing overlap upon inte- 
gration is \qn ~ qm\(- < 27r, where g is the local wavenum- 
ber within this domain. The probability that g„ would 
coincide with is 1/(A:£;£)''~^. The contributions of the 
non-zero overlaps add with random signs hence 



1/2 



X (*2^'^)vj, (C6) 



Assuming ergodicity w l/L'^^ and we get the same 
estimate as in the semi-classical procedure. 



Appendix D: The log-box ensemble 

The cumulative distribution function that corresponds 
to Eq. dMl) is 

In a; — In 



Prob(A < x) 



The algebraic, geometric and harmonic averages are 

xi - xo 



(Dl) 



In [xi/xq) 

{{x))g = ^XiXq 

{{x))h = ln(xi/a;o) 



XiXo 



(D2) 

(D3) 
(D4) 



Xi - Xo 

Note that for this distribution the median equals the ge- 
ometric average. The sparsity parameters are 

s = 2p 



P 



+ 1 

In p + In ( 1 



2p sinh 



(D5) 
(D6) 

(D7) 



where we defined p — (ln(a::i/a;o)) ^. If the distribution 
is very stretched reasonable approximations are 

s « 2p (D8) 

p —plnp (D9) 
For the VRH calculation 



Xl 



Xl 



p 



■ exp 



PQeUJ 



(DIO) 



For a rectangular F{uj) the VRH optimization is trivial 
and gives lu « x^^^ , leading to 



1 

?SLRT ~ - exp 
P 



1 



PQEi^C 



(Dll) 



For an exponential F(uj) the VRH optimization gives 



_ 1 
SsLRT ~ ~ exp 

p 



1 



1/2- 



(D12) 



^PQe^c, 

which is the same as in the traditional VRH optimization. 

Appendix E: The log-normal ensemble 

The cumulative distribution function that corresponds 
to Eq. (gni) is 



Prob(A < a;) = ^ + ^crf 



ln(x) — /i 



cr%/2 



The algebraic, geometric and harmonic averages are 

{{X))a = e^+-V2 

The sparsity parameters are 

s = 



The VRH estimate is 
x^ = exp 



+ aV2 erfinv ( 1 



(El) 

(E2) 
(E3) 
(E4) 

(E5) 
(E6) 

(E7) 
(E8) 



((a;))gexp 



llni'-] 6n^-lnln^'"' 



27r 



27r 



For a rectangular F{uj) the VRH optimization is trivial 
and gives uj « Xuj^ , leading to 



ffsLRT ~ q exp 



2V-lngln(fet^c)l (E9) 



For an exponential F{lu) the VRH optimization gives 



.9sLRT ~ q exp 



q exp 



— In q In 



-gluj'^lnq / V-4vr Inq 



^/-2\nq\n{gEU}c) 



(ElO) 



Due to the minority dominance the functional form is 
more robust compared with the log-box case. 
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FIG. 1: The driving induces transitions between levels E„ 
of a closed system, leading to diffusion in energy space and, 
hence, an associated heating. The diffusion coefficient De can 
be calculated using a resistor network analogy. Connected 
sequences of transitions are essential in order to have a non- 
vanishing result, as in the theory of percolation. 
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FIG. 2: Upper panel: The ergodicity of the eigenstates 
is characterized by the participation number (PN) which 
is calculated in various representations (see App. |B]). The 
bandwidth b of Vnm constitutes another measure for mix- 
ing. The clean, ballistic, diffusive and localization regimes 
(see Sec. Ill A[) are separated by vertical lines. Middle panel: 
The sparsity parameters (5, p and s) that characterize the 
perturbation matrix Vnm are plotted versus the disorder W . 
Lower panel: The scaled conductance in arbitrary units equals 
{(l^^mnP)). The Drude, the LRT and the SLRT results are 
displayed versus the strength of the disorder W . Inset: The 
same plot in the logarithmic scale. We see that in the ballistic 
regime the SLRT conductance becomes worse as the disorder 
becomes weaker, in opposition with the Drude expectation. 
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* resistors network calculation 
— VRH estimate 
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* resistor network calculation 
— VRH estimate 
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* resistors network calculation 
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FIG. 3: The SLRT suppression factor (?slrt versus the spar- 
sity parameter (p or q) for bi-modal distribution with rect- 
angular power spectrum (upper panel); (log-box distribution 
with exponential power spectrum (middle panel); log-normal 
distribution with rectangular power spectrum (lower panel). 
The resistor network and the VRH calculation were done 
for 100 realizations of 256 x 256 matrices with h — 10. In 
the big-modal case VRH is not satisfactory. In the log-normal 
case the VRH result is contrasted with the naive median based 
estimate. 
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FIG. 4: Upper panel: Time scales versus the disorder strength 
(see Sec llI Ap . Lower panel: Schematic illustration that sketch 
the dependence of the DC conductance on the strength of the 
disorder. It should be regarded as a caricature of Fig. [21 




FIG. 5: 

trix \Vnm 

regimes. The dashed lines correspond to the ballistic band 
width M — 10, which is associated with the time scale tL- 



Images of 100 x 100 pieces of the perturbation ma- 
P in the clean (W = 0.001) and ballistic (W = 0.1) 
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FIG. 6: The mean of the in-band values of the |wnmp ma- 
trix elements as a function of {n+m)/2 for different values of 
disorder. The pronounced modulation in the ballistic regime 
is an indication for texture. 
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FIG. 7: The ratio between the typical value [l^nm and the 
average value (|wnmp)cj that enter into Eq. ((6| and Eq. Q, 
respectively. In this plot the typical value is the median. The 
horizontal axis is n—m corresponding to the scaled frequency 
Qe(^. Note that both [|t;nmp]u and {\vnm\'^)Lj when plotted as 
a function of cj have a Lorentzian line shape. 
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(Anderson localization regime) 




FIG. 8: Schematic plot that illustrates the dependence of the 
average and the typical values of the matrix elements on the 
energy separation uj. These are labeled as 'LRT' and 'SLRT' 
respectively and compared with the semiclassical ('Drude') 
expectation. The plot refers to the Anderson regime. For 
a corresponding numerical illustration in the ballistic regime 
see Figdl 



